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Abstract. We present a novel method to recontruct the mass distribu- 
tion of galaxy clusters from their gravitational lens effect on background 
G^ I galaxies. The method is based on a least-x^ fit of the two-dimensional 

gravitational cluster potential. The method combines information from 
^ ■ shear and magnification by the cluster lens and is designed to easily 

incorporate possible additional information. We describe the technique 
and demonstrate its feasibility with simulated data. Both the cluster 
morphology and the total cluster mass are well reproduced. 
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G]N ■ 1. Introduction 
^ ■ 

Ph| It is of considerable interest for cosmology and the theory of structure formation in the 
O ' Universe to learn about the amount and the spatial distribution of mass in galaxy clusters. 

' While other methods to determine cluster masses depend on restrictive assumptions about 
ci • the symmetry and the equilibrium of the clusters, gravitational lensing is sensitive to the 
entire gravitating mass of the lens regardless of its composition and its dynamical state. 
I The drawback of lensing-based methods is that they determine the projected mass, thus 
' giving rise to ambiguities from projection effects, but these should be compared with the 
^ uncertainties of X-ray and dynamical mass estimates. 

Attempts to determine cluster mass distributions from lensing date back to the sugges- 
tion by Webster (1985) that clusters might act as efficient lenses on background populations 
of extended objects, and to the detection of a coherent shear pattern in the field of the clus- 
ter Abell 1689 by Tyson, Valdes, & Wenk (1990). Kochanek (1990) and Miralda-Escude 
(1991) discussed how parameterized cluster mass distributions could be constrained from 
observations of weak lensing in cluster fields. A technique for parameter-free cluster inver- 
sion was pioneered by Kaiser & Squires (1993), who found that the surface mass density 
of galaxy clusters can be derived in terms of a convolution of the observed shear pattern 
with a kernel describing the shear of a point mass. The feasibility of this and related 
approaches was then demonstrated by a number of authors (e.g. Bonnet, Mellier, & Fort 
1994; Fahlman et al. 1994; Small et al. 1995; Tyson & Fischer 1995), and Kaiser, Squires, 
& Broadhurst (1995a) and Bonnet & Mellier (1995) described techniques to accurately 
determine galaxy shapes from observed cluster fields. 
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The cluster inversion technique as proposed by Kaiser & Squires (1993) is designed 
for weak cluster lenses, and it causes unwanted boundary effects in small or irregularly 
shaped cluster fields. Schneider & Seitz (1995), Seitz & Schneider (1995a), and Kaiser 
(1995) extended the method into the nonlinear regime, and Schneider (1995), Kaiser et al. 
(1995b), Bartelmann (1995), and Seitz & Schneider (1996) described how the boundary 
effects can be removed. 

All cluster reconstruction techniques based on image distortions alone can determine 
the cluster mass distribution only up to a one-parameter family of linear transformations. 
This so-called mass-sheet degeneracy was originally found by Gorenstein, Falco, & Shapiro 
(1988) and discussed in the context of cluster lensing in the weak limit by Kaiser & Squires 
(1993) and generally by Schneider & Seitz (1995). The degeneracy can be broken if in- 
formation on the image magnification is employed, as suggested by Broadhurst, Taylor, &; 
Peacock (1995) and Bartelmann & Narayan (1995). 

Unambiguous determinations of cluster masses from weak lensing require the redshift 
distribution of background sources to be known. Apart from direct spectroscopy, which 
is hampered by the faintness of the background sources, lensing itself provides means to 
find faint-galaxy redshifts. Kneib et al. (1994) estimated arclet redshifts in the field of the 
cluster Abell 370, and Small, Ellis, & Fitchett (1994) discussed how faint-galaxy redshifts 
can be constrained from weak lensing by a number of clusters. Bartelmann & Narayan 
(1995) proposed an algorithm to simultaneously derive cluster mass distributions and the 
redshift distribution of faint galaxies. 

We propose here a novel cluster-reconstruction technique which is local, thus avoiding 
boundary effects by design, and which combines information on galaxy distortions and 
magnifications, thus breaking the mass-sheet degeneracy. The key idea is to reconstruct 
the two-dimensional Newtonian potential of the cluster lenses with a least-x^ approach. 
The least-x^ method lends itself to straightforwardly including additional information into 
the cluster reconstruction process. We aim at the gravitational potential rather than at the 
surface-mass density because it is the physical object underlying both lensing distortion 
and magnification, and because it should simplify comparisons of determinations of the 
cluster mass distribution from gravitational lensing with those from other means like, e.g., 
interpretations of the X-ray emission. 

This paper outlines the method and demonstrates its feasibility using simulated data. 
Extensions and a detailed discussion are postponed to a later, more technical paper (Seitz 
et al., in preparation). Section 2 describes the technique. In section 3, we present results 
from numerical simulations, and we summarize the paper in section 4. 



2. Outline of the method 

Gravitational lensing magnifies and distorts images of extended background sources. The 
local properties of the lens mapping are described by the Jacobian 




(1) 
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where x and y are dimensionless two-dimensional position vectors in the lens- and source 

planes, respectively. The convergence k is the scaled surface mass density of the lens, and 
7i 2 are the components of the shear. They are combinations of second derivatives of the 
(dimensionless) two-dimensional Newtonian potential of the lens, 

« = ^ (V',!! + "^,22) , 71 = ^ (V',!! - "^,22) , 72 = i^,12 ■ (2) 

For general reference on lensing, see Schneider, Ehlers, & Falco (1992) or Blandford & 
Narayan (1992). Image distortions, which can be quantified by, e.g., the quadrupole tensor 
of the surface-brightness distribution, measure the two-component quantity 

(Miralda-Escude 1991; Schneider 1995; Kaiser 1995). Since the ellipticity of an image is 
independent of its size, the ellipticity is unchanged if the Jacobian is multiplied by a factor 
A 7^ 0, which corresponds to transforming k and 7^ according to 

(1 - «) ^ A(l - «) , 7, ^ A7, . (4) 

The lens strength Qi above is manifestly invariant under this transformation, which reflects 
the mass-sheet degeneracy mentioned in the introduction. If the cluster has a critical curve, 
an ambiguity arises in the gi because of the parity change upon crossing the critical curve. 
An unambiguous measure of the ellipticity is then provided by the distortion 5^, 

2^.(1 -'^0 



(1 - + 7! + 7l 

(Schneider & Seitz 1995; Miralda-Escude 1991). 
The magnification jj, is given by 

^ = det-M = [(1 - ^)2 - ^2 _ ^2] -1 ^ (5) 

which scales with A~^ under the transformation (4). Measurements of the magnification 
can therefore be used to break the mass-sheet degeneracy, as Broadhurst et al. (1995) 
recognized. The magnification is accessible on a statistical basis by comparing the sizes of 
galaxies in cluster fields with those of galaxies of equal surface brightness in empty fields 
or by the change in number density of galaxies (Bartelmann & Narayan 1995; Broadhurst 
et al. 1995). 

We assume in the following that the sizes and shapes of galaxies in a cluster field 
and in an empty control field have been determined. The data are averaged over suitable 
regions of the cluster field to reduce the noise. These regions have to be small enough 
so that the properties of the lens can be assumed to be constant across each region, and 
large enough such that they contain a sufficient number of galaxies. Since the number 
of background galaxies is of the order of 10^ per square degree at i? ~ 27 (e.g. Tyson 
1994), this requirement is easily met. We do the smoothing with a Gaussian filter. The 
smoothing length is spatially constant and is varied later to minimize the obtained in 
the fit. 



4 



Maximum-Likelihood Cluster Recontruction 



For simplicity, we assume in the following that we already have the smoothed data on 
a regular grid covering the cluster field, i.e., that we have obtained smoothed measurements 
gi{k,l) and n{kj) of the two components of the lens strength and the magnification in 
each cell (/c, /) of the grid. For simplicity of notation, we use the inverse magnification 
r(/c, /) = ii~^{k, I) rather than the magnification. 

Of course, local estimates of the surface-mass density can be obtained directly once 
g and r are given at each grid cell (Broadhurst 1996). However, such an approach has 
the major disadvantage that it neglects the fact that convergence and shear are intrin- 
sically related fields. This relation provides important additional information, and thus 
considerably reduces the noise of the reconstruction. In our approach, the one physical 
field ip underlying all local lensing effects is reconstructed, and it is determined such as to 
optimize the global agreement of t/j with the data. 

We now want to determine least-x^ fits to the values of the potential ipik, I) in the cells 
{k, I) of the data grid. To do so, we replace the second partial derivatives of the potential 
in eq. [2] by their second-order finite-difference approximations. The resulting expressions 
can then be used to form estimators gi{k, I) and r{k, I) in terms of i/j for the lens strength 
gi{k,l) and the inverse magnification r{k,l) at the grid cell {k,l). The appropriate 
function then reads 

= E { ^2(^ ^9^{k, I) - Uk, l)f + [r(fc, /) - r(/c, , (6) 

where summation over the two components of gi is implied. Of course, depends on 
the values ip{m,n) of the potential at the grid points {m,n). The variances ag{k,l) and 
ar{k, I) of g and r can be estimated directly from the data. Since the derivatives 

d'ilj{m, n) 

are straightforwardly determined, can then be minimized numerically using a conjugate- 
gradient algorithm (e.g. Press et al. 1992, section 10.6). 

It is implicitly assumed in eq. [6] that g and r can be determined independently, which 
is the case in our simulation. In practice, however, g and r can be correlated especially 
for faint images. In that case, one would determine the local averages of g and r by giving 
lesser weight to fainter images, and modify the definition of such as to account for a 
possible correlation. 

The grid for the potential coincides with the grid of the smoothed data within the 
observed field, but it has to be extended by one further row or column of grid cells along 
each side of the field to allow the partial derivatives to be determined at any grid point 
where data arc given. We start the potential fit with il;{m,n) = for all {m,n). Since 
only second derivatives of the potential enter into the algorithm, we are free to add an 
arbitrary constant and a spatially constant gradient to the potential. We choose to fix the 
potential to ^ = at three corners of the field. The result is a least-x^ fit of the potential 
values i(^{m, n) on the grid, from which k and thus the surface-mass density can be derived 
via eq. [2]. In short, the x^ minimization determines the potential such as to optimize the 
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agreement between the measured ellipticities and sizes of the lensed galaxies and the shear 
and the convergence expected from the potential. 



3. Numerical simulations 

For testing the method, we use a numerically simulated galaxy cluster from the sample 
described in detail by Bartelmann, Steinmetz, & Weiss (1995). Briefly, the cluster was 
simulated within a COBE-normalized CDM cosmological simulation, taking the tidal forces 
of the surrounding structures into account. The cluster model we select is at redshift 
Zd = 0.16. It is resolved into ~ 17.000 particles whose line-of-sight velocity dispersion is 
~ 1200 kms-^ 

Background galaxies are assumed to be at a constant redshift of = 1. They are 
assigned random intrinsic ellipticity components (£i,£2) drawn from the probability dis- 
tribution 

TTaj [1 - exp (-l/c7|)] 
where jep = + el- For elliptical images, the modulus of e is given by 

II ^ ~^ 



a + h ' 



where a and h are the major and minor axes of the ellipse, respectively. We choose 
(Te = 0.15 (cf. Miralda-Escude 1991; Tyson & Seitzer 1988; Brainerd, Blandford, & Small 
1995). Further galaxy properties are chosen according to the galaxy model motivated 
and described by Bartelmann & Narayan (1995). They include luminosity and surface 
brightness, which allow to assign an intrinsic size to the sources. The logarithmic variance 
of the intrinsic galaxy radii is Z\ln(i?) a; 0.5, in good agreement with if^T measurements 
(Kaiser 1994, private communication). 

We simulate the lensing effect of the cluster on the sources, assuming a source density 
of 70 galaxies per square arc minute. In addition, we simulate an independent empty 
galaxy field to calibrate the intrinsic sizes of the sources needed for our algorithm. The 
field side length is 5', and the field is covered by a grid of 10 x 10 cells (or 12 x 12 cells for 
the potential). The simulated galaxy data in the lensed and the unlensed fields are then 
analyzed with the least-x^ potential reconstruction method described previously. We vary 
the smoothing length such as to minimize the obtained by the potential fit. Various 
results are presented in figures 1 and 2. The P^^ degree of freedom in this case is 1.08, 
with about equal contributions from the distortion and the size information. 

Figure 1 shows contour plots of the original cluster model in panel (a), the recon- 
struction in panel (b), the difference between these two in panel (c), and the potential 
in panel (d). A comparison between the upper panels shows that the cluster is very well 
reproduced except that the involved smoothing of the data tends to broaden the central 
mass peak. The agreement is quantified in panel (c), which shows that the residual is 
small and distributed fairly homogeneously across the field. In particular, there are no 
systematic deviations towards the field boundary. 



6 Maximum- Likelihood Cluster Recontruction 




Figure 1.- Four contour plots showing the original cluster model in panel (a), the reconstruction in panel 
(b), the difference between the two in panel (c), and the dimensionless two-dimensional potential in panel 
(d). Contours in panels (a) and (b) are spaced by 0.1 and the heavy contour follows k = 0.5. In panel (c), 
contours are spaced by 0.05 and the heavy contour follows Ak — 0. The potential is kept fixed at = 
at three corners. The heavy line in panel (d) follows the arbitrary contour ip = —5, and the contours are 
spaced by 1.5. The side length of the fields is 5'. 



Figure 2 displays two curves. The solid curve shows the total mass within radius 9 
from the center of the cluster relative to the input mass, i.e., the cumulative mass fraction 
within circles as a function of the circle radius. The dotted curve shows the mass fraction 
in annuli with outer radius 6 and width 0.'25. The solid curve starts at ~ 0.92, rises to 
~ 1.07, and falls towards ~ 1 at the field boundary. This shows that the total cluster mass 
enclosed by circles is reproduced to ~ 7% at any radius, and to higher accuracy at the field 
boundary. The mass fraction in annuli shows a somewhat larger fiuctuation of ~ ±10%. 
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Figure 2.- The reconstructed mass fraction within circles centered on the cluster center. The solid line 
shows the fraction of the cumulative mass within radius 0, the dotted curve the mass fraction within 
annuli. 9 is given in arc minutes. 

Both curves underestimate the mass in the cluster center, because the smoothing involved 

in the preparation of the input data tends to broaden the central mass peak, thereby 
shifting mass to larger radii. This local overestimate at intermediate radii is compensated 
by a slight underestimate at larger radii such as to reproduce the total mass, which is fixed 
by the overall magnification efi^ect of the cluster. 

4. Summary 

We have suggested a novel method to reconstruct cluster mass distributions from their 
gravitational lens effect on a population of background sources. Both distortion and mag- 
nification effects are included. The method is based on a least-x^ fit of the two-dimensional 
Newtonian potential of the lensing cluster. Reconstructing the potential is advantageous 
because it is the physical quantity underlying both the distortion and the magnification. 
The least-x^ approach has the advantages that it provides a local reconstruction technique, 
which renders it insensitive to boundary effects, and that it can be almost arbitrarily ex- 
tended to encorporate additional information. The combination of shear and magnification 
effects breaks the mass-sheet degeneracy inherent in all reconstruction algorithms based 
on ellipticity information alone. 
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Application of our method requires to measure not only the ellipticities of background 
galaxies, but also their intrinsic and their magnified sizes. Although the intrinsic dispersion 
of the sizes is large, the large number density of background galaxies supplies a sufficiently 
large data base to statistically extract reliable information on the magnification. The 
intrinsic sizes of the background galaxies are to be obtained from observations in empty 
fields. 

We have demonstrated using numerical simulations that the method accurately re- 
produces the shape and the total mass in the lensing cluster. The lensing cluster was 
numerically simulated within a CDM scenario, taking the tidal field of the surrounding 
matter distribution into account. The numerically simulated population of background 
galaxies was adapted to observations. It is designed to reproduce the observed intrinsic 
scatter in galaxy sizes and ellipticities and the number density of these objects. 

Various modifications and refinements are possible and can be included into the pro- 
posed technique. Among them are the effect of the lensing magnification on the number 
density of background galaxies (originally discussed by Broadhurst et al. 1995), and regu- 
larized maximum-likelihood techniques to avoid the (to some degree arbitrary) smoothing 
of the data and to ensure that the fitted data points are independent. Such modifications 
and extensions will be discussed in a following paper (Seitz, Schneider, Bartelmann, & 
Narayan, in preparation). 
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